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Low-lying A baryons with spin 1/2 are analyzed in full (unquenched) lattice QCD. We construct 
2x2 cross correlators from flavor SU(3) "octet" and "singlet" baryon operators, and diagonalize 
them so as to extract information of two low-lying states for each parity. The two-flavor CP-PACS 
gauge configurations are used, which are generated in the renormalization-group improved gauge 
action and the 0(a)-improved quark action. Three different /3's, /3 = 1.80, 1.95 and 2.10, are 
employed, whose corresponding lattice spacings are a — 0.2150, 0.1555 and 0.1076 fm. For each 
cutoff, we use four hopping parameters, (« V al, Ksea), which correspond to the pion masses ranging 
about from 500 MeV to 1.1 GeV. Results indicate that there are two negative-parity A states nearly 
degenerate at around 1.6 GeV, while no state as low as A(1405) is observed. By decomposing 
the flavor components of each state, we find that the lowest (lst-excited) negative-parity state is 
dominated by flavor-singlet (flavor-octet) component. We also discuss meson-baryon components of 
each state, which has drawn considerable attention in the context of multi-quark pictures of A(1405). 

PACS numbers: 
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I. INTRODUCTION 

A(1405) is one of the most interesting hadrons and at- 
tracting much interest from several view points. A(1405) 
is the lightest negative-parity baryon, even though it has 
one valence strange quark, which is heavier than the 
up and down quarks. Among the J p = l/2~ baryons, 
A(1405) is isolated from the others, much lighter than the 
non-strange counterpart A r (1535). It has no spin-orbit 
partner in the vicinity, as the lowest spin 3/2~ state is 
A(1520). Furthermore, the structure of A(1405) remains 
mysterious. On one hand, A(1405) is interpreted as 
a flavor-SU(3)-singlet three-quark state in conventional 
quark models. On the other hand, A(1405) could be in- 
terpreted as an antikaon-nucleon KN molecular bound 
state (B.E. ~ 30 MeV). The binding energy of KN im- 
plies a strong attractive force between K and N [HQ, 
which may cause a new type of dense hadronic matter, 
kaonic nuclei or kaonic nuclear matter [3 5]. We also 
expect that such KN bound states with large binding 
energies can be regarded as compact 5-quark states. The 
5-quark picture of A(1405) has advantages that all five 
quarks can be placed in the lowest-lying L = state to 
form a negative-parity baryon, and also that it requires 
no spin-orbit partner of A(1405). 

The property of A(1405) can therefore be an impor- 
tant clue to new paradigm in hadron physics. We here 
study properties of A(1405) using the lattice QCD for- 
mulation. Lattice QCD is nowadays employed as a very 
powerful tool for nonpcrturbative analysis directly based 
on QCD, and expected to cast light on the nature of 
A(1405) in a model-independent way. Though several 
lattice QCD studies on A(1405) have been performed so 
far p-|9|, most of them are based on quenched QCD 
and few lattice QCD studies succeeded in reproducing 
the mass of A(1405). Moreover, little has been discussed 
on the lattice about the possible mixing of flavor-SU(3)- 



octet and -singlet components induced by the symmetry 
breaking. 

Several possible reasons for the failure of reproducing 
A(1405) in lattice QCD were suggested through these 
studies, such as missing meson-baryon components due 
to quenching, exotic (non-3-quark type) structure of 
A(1405), or insufficiency of the lattice volume in the sim- 
ulations. Resolving such difficulties requires unquenched 
lattice QCD calculation on a larger lattice volume with 
varieties of interpolating operators. 

In this paper, we aim at clarifying the properties of 
A(1405) with two-flavor full lattice QCD, adopting the 
"octet" and "singlet" baryon operators to construct cor- 
relation matrices, which enables us to extract the low- 
lying spectrum as well as the mixing between octet and 
singles components in A(1405). 



II. LATTICE QCD SETUPS 

A. simulation conditions 

We adopt the renormalization-group improved gauge 
action and the O(a)-improved quark action. Simulation 
parameters are listed in TableQ] We adopt three different 
/3's, (3 = 1.80, 1.95 and 2.10, and corresponding lattice 
spacings are a = 0.2150, 0.1555 and 0.1076 fm [lj|, which 
arc determined so that the empirical p-meson mass is re- 
produced. The hopping parameters for strange quark k s 
are set to be 0.1431, 0.1393, and 0.1373 at /?=1.80, 1.95, 
and 2.10 so that the Kaon mass is reproduced. We note 
here that, if we adopt c/>-meson mass as an input, k s would 
be slightly different from those determined by the Kaon 
mass. The dynamics of A resonances would be closely 
connected to Kaon dynamics, and hence we adopt Kaon 
mass as an input in this study. We employ four different 
hopping parameters (K va i, K sea ) for each cutoff. Corre- 
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sponding pion masses range approximately from 500McV 
to 1.1 GeV at each (3. 



p 


csw 




x N t 


a [fm] 


N s a [fm] 




1.80 


1.60 


12 3 


x 24 


0.2150 


2.580 


0.1431 


1.95 


1.53 


16 3 


x 32 


0.1555 


2.488 


0.1393 


2.10 


1.47 


24 3 


x 48 


0.1076 


2.582 


0.1373 



TABLE I: List of the present simulation parameters. /3 and 
csw are the parameters for the renormalization-group im- 
proved gauge action and the 0(a)-improved quark action. N s 
and Nt denote the spatial and temporal lattice extents, and 
the lattice spacings a and the strange-quark hopping param- 
eters k s are determined so that the empirical p-meson mass 
and Kaon mass are reproduced • 



B. baryonic operators 

In order to extract the low- lying states in S = — 1 and 
isosinglet channel, we construct 2x2 cross correlators 
from the following "singlet" and "octet" operators, 

m{x) = Ai(s) = -±=e abc {u a (x)[d Tb (x)C l5 s c {x)\ 
+d a {x)[s Tb {x)C l5 u c {x)}+s a (x)[u Tb {x)C l5 d c {x)}} (1) 

m (x) = A 8 (x) = ±=e abc {u a (x)[d Tb (x)C l5 s c (x)} 

+d a {x)[s Tb (x)C l5 u c (x)} ~ 2s a (x)[u Tb (x)C l5 d c (x)}} (2) 

It is easy to check that 771(2;) (r)s(x)) belongs to the sin- 
glet (octet) irreducible representation of the flavor SU(3). 
The cross correlator of 771 and r/s vanishes in the fiavor- 
SU(3) symmetric limit (m u — nid = m s ). We define 
cross correlators as 

M(x,y) IJ = Tr^± ( m (x)fjj(y)), (3) 

where the trace is taken over spinor indices. The parity- 
projection operator, P(±) = "^p-, is inserted to separate 
the positive and negative parity states. In the actual sim- 
ulation, we adopt point-type operators for the sink, r](x), 
and extended operators, which are smeared in a gauge- 
invariant manner, for the source, fj(jj = 0). Smearing 
parameters are chosen so that root-mean-square radius 
is approximately 0.5 fm. 



C. diagonalization of cross correlators 

We here consider a general situation where we have 
a set of TV independent operators, 77f nk for sinks and 
?7j 10 ^ for sources in order to construct correlation matri- 
ces M.u(t) = (?7f nk (i)?7j rC ^(0))j which can be decomposed 



into the sum over the energy eigenstates \i) as 
Mu(t) = M IJ (t,0) = (r 1 r*(t)r ] r f m 

= E( C S t nk)/^(t) lJ (C src ) jJ 

i,j 

= (Cl k A(t)C SIC ) u , (4) 

where the small letters (ij) are the indices for the inter- 
mediate mass eigenstates and 

(CU)n = (™\vT k \i), (C TO ) 3 -/ = (il^rVac), (5) 

are the coupling matrices. A(t) is a diagonal matrix given 
by the energy eigenvalue, Ei, of the state i. 

A(% = 5^^. (6) 

From the product 

M-\t + l)M(t) = C£A(-l)C m , (7) 

we extract the effective eigen-energies (or the effective 
masses after the zero-momentum projection) {Ei} at the 
time slice t as the logarithm of eigenvalues {e Ei } of the 
matrix M^(t + l)M(t). 

Besides overall constants, (C src ) _1 and (C^jJ -1 can be 
obtained as right and left eigenvectors of M~ 1 (t+l)A4(t) 
and Ai(t)Ai(t + l) -1 respectively, since 

M-\t + lVW^XCWc)- 1 = (C src )- 1 A(-1) (8) 

and 

(Cl^M^Mit + l)- 1 =A(-l)(C| nk )- 1 (9) 

hold. One can also directly construct optimal source and 
sink operators, O s * c ^ and Gf lk , which couple dominantly 
(solely in the ideal case) to i-th lowest state, as 

0-t = ^ r t (Csrc) -i (10) 
j 

and 

When an optimal correlator which corresponds to a 
single-state propagation is needed, one can extract it by 
sandwiching M(t) between (Cj nk ) _1 and (C^)" 1 since 

(C s t nk )- 1 M(t)(C src )- 1 =A(t) (12) 

is diagonal. We note here that, if the correlation ma- 
trix A4(t) is hermitian, one can determine (C src ) _1 and 
(Csnk) -1 up to overall phase factors so that Eq. (fl"2|) is 
satisfied. 

These prescriptions are valid only in the i-range where 
the number of the relevant physical states is no more than 
N, since we prepare a set of just ./V independent opera- 
tors. In practice, higher excited states may come in, and 
such contamination manifests itself as a i-dependence 
in several "constants" mentioned above. In order to 
avoid the contamination of the excited states, we iden- 
tify a t-window in which the "constants" are (almost) 
t-independent. All the physical quantities are extracted 
within such a i-window. 
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D. different choices in formulation 

We here discuss the difference in formulation. Ingen- 
cral, we solve generalized eigenvalue problems [ill fl2l]. 

M-\t')M{t) = {C~l)h(t - t')(C SIC ), (13) 
M{t)M-\t') = (CL)A(i - t')(cUr\ (14) 

where we have ambiguity in the choice of t' . As long as 
a system is dominated by the states as many as the rank 
of correlation matrices, choice dependence vanishes. In 
practice, we may have excited-state contaminations, and 
choice dependence can emerge. In our case, \t— 1'\ is fixed 
to 1, and we take as large t as possible: We seek plateau 
regions and fit data in that window. The authors in 
Ref. [l3[ investigated possible subleading contaminations 
by higher excited states and the dependence on the choice 
of t and t' . 

To estimate the systematic errors coming from such 
ambiguity, we try four different choices, \t — t'\=l, 2, 3, 
and 4, and see the dependence. Namely, we solve the 
eigenvalue problems 

M(t)M-\t + m) = (C s t nk )A(-m)(C s t nk )- 1 , (15) 

with m < 4. We plot in Fig. Q] (C snk ) i/(C snk ) 8, 
which can be a measure of flavor contents in an eigen- 
state and will be discussed later. They are obtained at 
K SC a,K va i = 0.1367 at /3 = 2.10. Though they are all t- 
independent in an ideal case, there remain considerable 
contaminations especially with small m. On the other 
hand, it reaches plateau region much faster with larger 
m. However larger m results in larger statistical fluctu- 
ations and diagonalization sometimes fails at relatively 
small t. In our case, as long as we identify plateaus prop- 
erly, resulting values are all consistent with each other. 

III. LATTICE QCD RESULTS 

A. hadronic masses 

In Fig. [2j as typical results, effective masses are plot- 
ted for the ground and the lst-excited states obtained 
from positive- and negative-parity projected data with 
K S oa,K va i = 0.1367 at j3 = 2.10. We can find i-windows 
in which the effective masses exhibit plateaus. We ex- 
tract eigen-energics by fitting lattice QCD data in the 
t-range where both of the ground- and lst-excited states 
exhibit plateaus. In Table [Til we list all the hadronic 
masses extracted in this study. Figs. [3] and [4] show the 
eigen-energies in the positive- and negative-parity chan- 
nels, respectively, plotted as functions of the squared pion 
mass m%. Filled circles denote the energies of the ground 
states, and open squares those of lst-excited states. The 
solid curves represent quadratic fits as a function of 
squared pion mass mj. Two solid lines at the verti- 
cal axes indicate the masses of A(1115) and A(1600) in 
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FIG. 1: (C sn k)oi/(C' an k)o8 obtained by solving eigenvalue 
problem, M{t)MT^{t + m) = (C s t nk )A(-m)(C s t nk )- 1 with 
m < 4 are plotted as a function of t. They are obtained 
at fi; SC a, ftvai = 0.1367 at ft = 2.10. The horizontal lines de- 
note the fitted value and the corresponding error evaluated 
by fitting plateaus. 



Fig. [3] (positive-parity states), and A(1405), A(1670) and 
A(1800) in Fig. |4] (negative-parity states). 

The obtained masses of the positive-parity ground 
state agree very well with the mass of the ground- 
state A(1115) at all three /3's. On the other hand, the 
lst-excited state in this channel lies much higher than 
A(1600), which is the 1st excited state experimentally 
observed so far. The same tendency was reported in 
Ref. [H, and the situation is similar to the case of the 
Roper resonance, which is the non-strange SU(3) part- 
ner of A(1600) [jj]. 

In the negative-parity channel (KigHJ) , the ground- and 
the lst-excited states always have close energies at all the 
/3's and k's. The eigen-energies have similar quark-mass 
dependences, and the mass splittings are almost quark 
mass independent. The chirally extrapolated values both 
lie around the mass of A(1670) rather than A(1405). Sim- 
ilarly to the previous studies, we do not reproduce the 
mass of A(1405) in our calculation. While in quenched 
simulations in Refs.0, 0] such failure was regarded as an 
evidence of possible meson-baryon molecule components 
in A(1405), our present unquenched simulation contains 
effects of dynamical quarks and thus should incorporate 
meson-baryon molecular states. We will come back to 
this point later in section V. 



B. lattice-spacing dependences 

Fig. [5] shows the chirally extrapolated masses of the 
lowest two states in each parity v.s. the lattice spacing 
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FIG. 2: The effective mass plots for ground and lst-excited 
states in the positive- (upper panel) and negative-parity 
(lower panel) channels, which are obtained with 
0.1367 at ft = 2.10. The horizontal lines are fits of the eigen- 
energies in the plateau region for each state. 



a. Except for the lst-excited positive-parity state at the 
largest a, all the data show little lattice-cutoff depen- 
dence. The horizontal lines indicate the masses of the 
experimentally observed A states. The positive-parity 
ground state agrees well with the mass of A(1115), while 
the lowest two negative-parity states lie around A(1670). 
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C. possible contaminations by scattering states 

Since our calculations contain dynamical u and d 
quarks, scattering states could come into spectra. NK 
and S7T thresholds are drawn in Fig. [3] and Fig|4] as 
dashed lines. The energy of the lst-excited state at 
the lightest u- and d-quark masses could be contami- 
nated by such scattering states. Several strategies have 
been adopted to discriminate resonances from scattering 
states: study of volume dependence of the masses, and 
the spectral weights or boundary condition dependence 
of the eigen-energies. They were tested and applied to 
the pentaquark baryon p^l - l2fj| |. 

There are two meson-baryon channels relevant in the 
present calculation, 7rE and KN. In terms of the valence 
quark contents, we have 5 different thresholds (2 for KN, 
3 for 7rS). In this paper, in order to distinguish resonance 



FIG. 3: Masses of the lowest two positive-parity A states as 
functions of the squared pion mass. The filled circles (open 
squares) denote the masses of lowest (lst-excited) state. Two 
solid curves represent quadratic functions used in the chiral 
extrapolation. Two dashed lines indicate the 7rE and the KN 
thresholds in the presence of the relative momentum p = 
Two solid lines on the vertical axes show the experimentally 
observed masses of A(1115) and A(1600). 



states from all these scattering states, we impose the fol- 
lowing boundary condition on the quark fields, 

ip(x + L) = e% ni ip(x). (16) 

Under such boundary condition for quark fields, a 
hadronic state 4>zk+n{x) which contains 3/c + n valence 
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(j.L. 


0.455(16) 
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Ci A 1 A / 1 O A 

0.919(18) 


2.10 0.1357 0.630(1) 
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0.1374 0.424(1) 
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0.565(1) 
0.491(1) 
0.435(1) 
0.363(1) 


1.194(05) 
1.016(05) 
0.902(05) 
0.734(06) 


1.131(06) 
0.990(05) 
0.907(05) 
0.784(05) 


C.L. 


0.312(9) 


0.579(32) 


0.674(26) 






Et 


Eo 


EI 


1.80 0.1409 2.201(03) 
0.1430 2.013(04) 
0.1445 1.836(04) 
0.1464 1.533(07) 


3.53(34) 
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2.10 0.1357 1.133(09) 
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1.526(15) 
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1.325(15) 


1.384(13) 
1.246(09) 
1.158(10) 
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1.437(12) 
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0.977(31) 
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TABLE II: Hadronic masses (in lattice units) obtained for 
four different k's at each /3. The m w , m,K, m,N, and ms denote 
the masses of pion, kaon, nucleon, and E baryon, respectively. 
Ef represents the eigen-energy of the i-th state in positive 
(+) or negative ( — ) parity channel with 1 = and S = — 1. 
"C.L." in each block shows the results extrapolated to the 
chiral limit with quadratic functions. 



FIG. 4: Masses of the lowest two negative-parity A states as 
functions of the squared pion mass. The filled circles (open 
squares) denote the masses of the lowest (lst-excited) state. 
Two solid curves represent quadratic functions used in chi- 
ral extrapolation. Two dashed lines indicate the 7rE and 
KN thresholds. Two solid lines on the vertical axes show 
the experimentally observed masses of A(1405), A(1670) and 
A(1800). 



quarks obeys 

fo k+n {x + L) = e§ n7r V3fc+n(z), (17) 

and can have spatial momenta, 

27r 2mr 
Plat = -j-m + --j- (m G Z) (18) 



As a result, only states which consist of 3k valence 
quarks can be zero-momentum states. Since other quark 
combinations inevitably have non-vanishing spatial mo- 
menta, their energies are raised up. As long as we employ 
three-quark operators for baryon creation/annihilation, 
hadronic states appearing in scattering states should con- 
tain one or two valence quark(s), since sea-quark pairs 
themselves cannot carry flavors. 

Now we turn back to our full-QCD calculations. We 
impose the following boundary condition to all the va- 
lence quark fields (q = u, d, s); 

q(x + L) = ei ni q(x). (19) 

We plot in Fig. |6] the eigen-energies under the peri- 
odic and the twisted boundary conditions, which are ob- 
tained at (3 = 2.10 and n = 0.1382. The open and filled 
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n 

P K 


9a x 10^ 


9i x !0 2 


9o x w * 


9i x 10 2 


1.80 0.1409 
0.1430 
0.1445 
0.1464 


4.91(303) 
0.12(013) 
-3.82(195) 
0.10(126) 


0.30(01) 
0.02(00) 
-0.33(01) 
-1.39(06) 


-10.7(036) 
-0.64(035) 
15.2(023) 
35.7(077) 


6.46(291) 
0.40(027) 
-10.5(018) 
-21.9(063) 




—8.36(800) 


— 1.66(27) 


65.3(116) 


AO 1 / l 1 0\ 

—43.1(113) 


1.95 0.1375 
0.1390 
0.1400 
0.1410 


7.22(066) 
0.58(029) 
-2.46(022) 
-7.02(081) 


0.46(01) 
0.11(00) 
-0.32(00) 
-1.14(01) 


-33.4(145) 
-2.99(141) 
11.4(013) 
26.3(031) 


22.2(109) 
4.25(167) 
-8.83(120) 
-15.2(025) 


CL. 


—9.12(323) 


— 1.81(24) 


39.9(049) 


— 19.8(072) 


Z.1U U.loO/ 

0.1367 
0.1374 
0.1382 


1U.U ( (128 J 

5.17(031) 
-0.77(010) 
-5.44(112) 


U. /8(Uo J 

0.37(01) 
-0.09(00) 
-1.11(04) 


— zz.o^Uoo J 
-17.9(061) 
2.19(051) 
19.8(036) 


12.4(^UoD J 

8.21(367) 
-1.26(044) 
-15.2(040) 


C.L. 


-14.4(052) 


-1.81(47) 


42.5(103) 


-33.3(017) 



TABLE III: The ratios gf (multiplied by 10 2 ) denned in 
Eqs. (|23|I - (|26|I . "C.L." shows the values extrapolated to the 
chiral limit using quadratic functions. 



circles (squares) are for the ground (lst-excited) states. 
The lower and upper solid lines respectively represent 
the threshold energies of E-7T and NK states with nor- 
mal/twisted boundary condition. The threshold energies 
are raised up in the case of twisted boundary condition, 
whereas the lowest two eigen-energies remain unchanged. 
Thus we conclude that the observed states are insensitive 
to boundary conditions, and contaminations from meson- 
baryon scattering states are small. Namely the results 
indicate that these states are regarded as compact reso- 
nance states. 

Similar results are obtained also for the positive-parity 
states. Their energies are insensitive to the boundary 
conditions. (The energies of the positive-parity lst- 
excited states at the smallest n are higher than p-wave 
thresholds (with relative momenta p = =£) of £7r or NK 



states.) 



IV. FLAVOR STRUCTURES 

The chiral unitary approach [2l| has suggested that 
A(1405) is not a single pole but a superposition of two 
independent resonance poles. The structure of A(1405) 
is now attracting much interest, and desired to be clari- 
fied in a model independent manner. We investigate the 
flavor structure of the ground and the lst-excited states 
obtained from the cross correlators of two operators. 

In order to clarify the flavor structure of the low-lying 
states, we evaluate the overlaps of the obtained states 
with the source and sink operators. We define the mag- 
nitudes of the singlet (1 = 1) and the octet (1 = 8) 
components in the i-th state (i = for the ground state 
and i = 1 for the 1st excited state) by 
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FIG. 5: Masses of the ground and the lst-excited states in 
the chiral limit for the positive (upper) and negative (lower) 
parity A states plotted as functions of the lattice spacing a. 
The horizontal lines indicate the masses of the experimentally 
observed A states. 



Actually, we can only obtain the ratio between 

(i|T7+|vac) = (CV (21) 

and 



(C) ia , 



(22) 



since the overall factors in C src cannot be determined in 
our setups. We evaluate go and g\ defined as 



.9n = Cos/Cqi 
9i = Cii/Cis, 



(23) 
(24) 



Both and g-j" vanish when the SU(3) f symmetry is ex- 
act, showing that the ground (1st excited) state is purely 
flavor singlet (octet) in the limit. For the positive-parity 
states, we similarly define 



So = Coi/Cos 



9l - 



(25) 
(26) 



(C). 



i- 



(20) 



In this case, we exchange the denominator and the nu- 
merator since the ground state is flavor octet and the 1st 
excited state is flavor singlet in the SU(3)f limit. 

We show gQ and gf 1 as a function of pion-mass squared 
for each lattice cutoff in Fig. [7] and Fig. [5] . As is ex- 
pected from the symmetry, such mixing coefficients cross 
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■ 1 -st state 
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FIG. 6: Lowest two eigen-energies in the (J P , S) = (1/2 - ,— 1) 
channel at /3 = 2.10 and k = 0.1382 under the normal (pe- 
riodic) and the twisted boundary conditions. The circles 
(squares) are for the ground (lst-excited) states. The lower 
and upper solid lines respectively represent the threshold en- 
ergies of E-7T and NK states evaluated with normal/twisted 
boundary condition. 



zero at the flavor symmetric limit (k u = = k s ). The 
data show smooth quark-mass dependence toward the 
chiral limit in all the lattice-cut-off cases and the depen- 
dences are almost lattice-cut-off independent. We hnd 
that the magnitude of operator mixing gets larger for 
larger flavor-symmetry breaking. It is interesting that 
the lst-excited state in positive-parity channel is domi- 
nated by the flavor-singlet component showing almost no 
contaminations of octet components. 

We perform quadratic fits of these mixing coefficients 

as 

go(ml) = a 2 fiml + ai,om» + a ,o (27) 
9i(ml) = a 2 . 1 ml + a ltl m^ + a ,i (28) 

and list the extrapolated values in Table IIIII Even in 
the chiral limit, gi (i = 1,2) in negative-parity channel 
remain less than 0.6, which indicates that the ground 
(lst-excited) state is dominated by the singlet (octet) 
component in the present quark-mass range. 

On the other hand, the mixings of the singlet and octet 
components are generally weaker in the positive parity 
A's. It is natural because the splitting of the two states 
are large in the SU(3)p. Nevertheless, it is interesting to 
see that the ground state, corresponding to A(1115), con- 
tains significant (about 10%) flavor singlet component in 
the chiral limit, which was not expected from the simple 
quark model picture. 

In practice, we need to interpret such overlap coeffi- 
cients into physical amplitudes. In the Euclidean-time 
range we consider 



r\ x | vac) 
77s | vac) 



co|0)+siCi|l) 
g o c o \0) +ci|l) 



(29) 
(30) 
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FIG. 7: gl are plotted as a function of m n . Solid lines denote 
quadratic functions used in chiral extrapolation. 



expanded in terms of the singlet and octet basis states as 



|0> 

ID 



MoilAO 
Mii|Ai) 



M 08 |A 8 ) 
M 18 |A 8 ) 



(31) 
(32) 



are satisfied. We here assume that each state can be 



We note that |A)'s in this expression are no longer en- 
ergy eigenstatcs in QCD, and they are restricted to el- 
ements that interpolating operators have overlaps with; 
| A) cx jy^vac). (Actually we can introduce other singlet 
and octet bases, for example, whose spatial distributions 
are different. If we adopt more interpolating fields such 
as those with spatial derivatives or with exotic spinor 
structures, we can introduce much more basis states in 
this argument.) 
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Using the normalization relations, 
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FIG. 8: g i are plotted as a function of m%. Solid lines denote 
quadratic functions used in chiral extrapolation. 



We here extract My's, using the conditions, 
(^8 hi | vac) = and (Ai|?7 8 |vac) = 0, giving 

M 08 



M 18 
Mi 
Mo 



-9i R, 
i = -g R- 1 



(33) 
(34) 



with R — ci/crj. R 2 can be estimated from the 2-point 
correlators, since optimized operators satisfy 

{m - ffi77s)|vac) = co(l - 5o5i)|0) (35) 
(»7s ~ SWi)|vac) = ci(l - 3o5i)|!>- (36) 



Mp 2 ! + M 2 S = 1 



J 08 



M 2 1 + M 2 S 



(37) 
(38) 



we determine M 2 - , the amplitude of j-plet component in 



i-th state as follows: 



M 2 ! 



M 2 m 



Ml 



i-gfR 2 

1 - (.90.91 ) 2 
9l{R 2 -9%) 
1 - (203i) 2 

go 2 (l-g^ 2 ) 
R 2 {l-(g gi) 2 } 

R 2 -9o 

R 2 {1 - (gogi) 2 } 



(39) 
(40) 
(41) 
(42) 



The amplitudes M 2 t obtained on the finest lattice are 
plotted in Fig. |H1 In the present quark-mass range, 
ground (lst-excited) state is dominated by flavor-singlet 
(octet) components. 
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y obtained at /3 = 2.1 are plotted as a function 



FIG. 9: Mfj 
of squared pion mass 
Mf 8 = l. 



Mij are normalized so that Mf\ + 



So far, we have argued the internal structures from 
the view point of flavor contents. It would be interest- 
ing and useful to interpret these flavor-mixing amplitudes 
in terms of a hadronic basis, since hadronic-basis ex- 
pansion of unstable particles could be directly compared 
with experiments and sometimes beneficial for model cal- 
culations. Each flavor element can be decomposed into 
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hadronic bases as 



|Ai) 
|A 8 ) 



W; 



NK.X 



w> 



NK.8 



NK) 
NK) 



(43) 
(44) 



In this paper we restrict ourselves only to the NK and E-7T 
bases. The SU(3) Clebsh-Gordan coefficients tell us [22| 



|Ai) cx 



(y/2\NK) + V3\Xn) - | At;) - V2\EK) 



V8 
|A 8 )cx 

3D 



(45) 



'20 
9f 



V2\NK) - yi2|E7r) - VZ\Arj) - V2\EK) 



V6\NK) + y/6\EK 



and the following expressions follows; 



(46) 



W 



NK.l 



W N k,S 



TF 



E7r,: 




Here, gu and <?i? are unknown coupling constants 
corresponding to two possible SU(3) constructions of 
the meson-baryon octet states, Tr ({.£?,£?} M) and 
Tr (\B, B] M), respectively. We finally reach the proba- 
bilities of NK and T,tt states in i-th state, P l NK and 



We focus on the meson-baryon components of the low- 
lying states in negative-parity channel, since they are of 
interest here. We show in Fig. [TUl PL~(a) and P|. 7r (a) 
normalized as P^^(a) + P^ n (a) = 1 as a function of 
a. Though g^s are determined from those obtained at 
f3 = 2.10, the behaviors do not change significantly even 
if we adopt extrapolated values on coarser lattices. Now 
that the only unknown parameter is a representing gu 
and gF in the J p = 1/2 - channel. Unfortunately, the 
couplings, gu and gp, should be determined up to in- 
cluding their signs in order to estimate the meson-baryon 
components in each state in this strategy. They can in 
principle be determined from lattice QCD computations 
and their determination is left for forthcoming studies. 
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FIG. 10: P^jf(a) and P^(a) of the ground state (upper 
panel) and of the 1st excited state (lower panel) as functions 
of q. They are determined from the g~'s obtained at /3 — 2.10 
(the finest lattice). 



Defining a by 

(cosa,sina) = ( 9D 9F ) (53) 

W9d+9f V9d+9fJ 

we obtain P l N ^(a) and P^{a) as functions of the angle 
a, where we substitute go and g± extrapolated to the chi- 
ral limit at (3 = 2.10 (the finest lattice). We simply adopt 
R 2 = 1, since R 2 does not deviate from unity, and is al- 
most quark-mass independent; R 2 =1.002(29), 1.015(21), 
1.030(29), and 1.024(29) at K u , d =0.1357, 0.1367, 0.1374, 
and 0.1382 at this f3. The relative signs in My are de- 
termined so that MiiMis(C)n(C)i2 > 0. 



V. CONCLUSION AND DISCUSSIONS 

An unquenched lattice QCD study for the low-lying A 
baryon states in the S = —1, 1 = and J = 1/2 channel 
has been performed, focusing especially on the masses 
and structures of the two lowest negative-parity A's. We 
have constructed 2x2 cross correlators of the "octet" and 
"singlet" baryonic operators, and diagonalized them so 
that we can properly extract the information of the low- 
lying two states. In our measurements, we have adopted 
the 2-flavor gauge configurations generated by CP-PACS 
collaboration with the renormalization-group improved 
gauge action and the C(a)-improved quark action. We 
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have performed simulations at three different /3's, /3 = 
1.80, 1.95 and 2.10, whose corresponding lattice spacings 
are a = 0.2150, 0.1555 and 0.1076 fm, employing four 
different hopping parameters (K V ai,K sea ) for each cutoff, 
corresponding to the pion masses ranging from 500 MeV 
to 1.1 GeV. 

It is shown that the mass of the ground-state A(1115) 
is reproduced well at all three /3's, while the lst-excited 
positive-parity A state lies much higher than the experi- 
mentally observed A(1600). The same tendency was re- 
ported in Ref. Q , and also in the case of the Roper reso- 
nance, iV*(1440), which is the non-strange SU(3) partner 
of A(1600) [3]. It should also be noted that no 7r£ or 
KN scattering states are seen in the calculation. 

Our results for negative-parity A states indicate that 
there are two l/2~ states nearly degenerate at around 
1.6— 1.7GeV, while no state as low as A(1405) is observed. 
We have revealed the flavor structures of these states 
from the lattice data for the first time. It is found that 
the lowest negative-parity A state is dominated by flavor- 
singlet component. The second state, which is less than 
100 MeV above the ground state, is predominantly flavor 
octet. Thus we find that the first two negative-parity A's 
have different flavor structures. 

They, however, come significantly heavier than the ex- 
perimentally observed lowest-mass state, A(1405). One 
naive possibility is that the obtained states do not cor- 
respond to the physical A(1405), but describe excited A 
states. In fact, the 2nd and 3rd negative-parity states 
lie at 1670 MeV and 1800 MeV, both of which are three- 
star states in the Particle Data Group classification. The 
present lattice data are consistent with these excited 
states. In the non-relativistic quark model approach, 
each of these states is classified as a flavor octet P-wave 
baryon. We, however, have shown that the lower state is 
dominated by a flavor-singlet component. So, the results 
here predict one flavor-singlet state and one flavor-octet 
state in the vicinity. The reason for missing A(1405) state 
will be its poor overlap with three-quark operators. In 
fact, the inclusion of dynamical quarks does not strongly 
enhance the signals of possible meson-baryon scattering 
states. (See the later discussions.) If A(1405) is predom- 
inantly a meson-baryon molecular state, such overlaps 
would be naturally small. 

The other possibility is, of course, that the lowest 
and the 2nd-lowest states describe physical A(1405) and 
A(1670) but the masses have been overestimated. Con- 
sidering that our simulation contains two-flavor dynami- 
cal quarks, the failure of obtaining a light A state could 
be attributed either to (1) strange-quark quenching, (2) 
insufficient lattice volume or (3) lack of chiral symmetry. 
In fact, (1) seems to make the masses of octet baryons 
in positive parity channel slightly (< 10%) overestimated 
in the present setups [To|. On the other hand, the defi- 
ciencies, (2) and (3), may cause the lowest state not prop- 
erly reproduced, supposing that the main component of 
A(1405) is a meson-baryon molecular state. In order to 
check whether this conjecture is correct, simulations with 



light dynamical quarks (m„- <C 500 MeV) and larger vol- 
ume (L 3> 2.5 fm) will be required. 

Upon the above conjecture, we can further consider one 
interesting scenario that these states both correspond to 
the physical A(1405). Then, the results may support 
its double pole structure proposed by the chiral unitary 
approach [21] . In our results, the lowest two states are 
almost degenerate at all the /3's (lattice spacing) and k's 
(quark masses). Namely, the obtained two states are the 
signature of the double-pole resonance, but the mass has 
not yet been reproduced because of the deficiency stated 
above. 

It is also important to note the missing meson-baryon 
scattering states. Despite that dynamical up and down 
quarks are included and the meson-baryon thresholds ap- 
pear around/below the obtained eigen-energies, we have 
found no clear signal of the meson-baryon scattering 
states. This fact is supported by the observation that 
the effective mass plots under the normal and twisted 
boundary conditions show no prominent differences. We 
show in Fig[TT]the effective mass plots of the ground state 
negative-parity A under the normal and twisted b.c.'s. 
They are obtained on the 24 3 x 48 lattice with the largest 
hopping parameter. Thus we conclude that no scattering 
states appear in the present spectrum. Also in Ref. [23| . 
in which excited-state nucleon spectrum was systemati- 
cally and extensively investigated with two-flavors of dy- 
namical quarks, no clear signal of scattering states was 
found and the importance of multi-quark operators was 
raised. 

A similar situation can be found in the computation of 
the Wilson loops, whose expectation values give us the 
potential between a (heavy fundamental) quark and an 
antiquark. In the presence of dynamical quarks, such an 
interquark potential should saturate and flatten at some 
interquark distance, where the confining string is broken 
and a quark- antiquark pair is created. However, no one 
has ever observed successfully such a "st ring breaking" 
effect in the Wilson loop computations [24[ [25| . One 
possible reason for this phenomenon is that the Wilson 
loop itself has poor overlaps with such broken strings. 

In the A spectrum, we may similarly conjecture that 
the scattering states or broad resonances, which are sen- 
sitive to boundary conditions, have little overlap with the 
3-quark interpolating field operators. More detailed in- 
vestigation would be needed for clarification of scattering 
states purely induced by dynamical quarks. 

The flavor structures of the A states have been very 
well clarified using the variational method. The octet and 
the singlet components are mixed when the flavor-SU(3) 
symmetry is broken. Actually, the ground (lst-excited) 
state is dominated by singlet (octet) component, and the 
contamination by another representation is at most 20% 
(5% when squared) in our present analysis. From these 
findings, we expect that the flavor-SU(3) symmetry is 
not largely broken. A similar conclusion was also derived 
for the study of the meson-baryon coupling constants in 
lattice QCD 
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FIG. 11: Effective masses of the negative-parity A ground 
state under the normal (periodic) and the twisted boundary 
conditions. 



Because the SU(3) breaking effect seems small, the 
analyses without SU(3) mixings adopted so far 
make sense to some extent. The mixings, however, get 
larger towards the chiral limit, variational analyses could 



be essentially needed when we adopt much lighter quarks. 
We also find the meson-baryon contents in each state 
strongly depend on the meson-baryon couplings and their 
signs. Precise determination of the meson-baryon con- 
tents will require reliable determination of the couplings 
up to signs, for which further lattice QCD calculations 
may be helpful. 
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